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Abstract 


Resistive  instabilities  with  azimuthal  mode  numbers  m  =  0  and  1  are 
studied  on  intense,  relativistic  electron  beams.  An  exact  cold-fluid  model 
of  the  beam  is  used,  and  the  background  is  modeled  by  a  scalar  conductiv¬ 
ity.  The  effects  of  the  beam  profile,  the  width  of  the  conductivity  and 
return  current  profiles,  and  the  magnitude  of  the  return  current  fraction 
are  considered  for  beams  in  the  10-100  kA  range  with  y  =  100.  Some  results 
for  modes  with  m  =  2  and  3  are  also  presented. 


I 


I 

Introducti on 

The  most  serious  instabilities  of  intense  electron  beams  propagating  in 
a  highly  collisional  plasma  appear  to  be  the  resistive  modes1-3  (juj  << 

4ito,  where  w  is  the  frequency  of  the  mode  and  o  is  the  background  conduc¬ 
tivity).  In  this  report,  we  pay  particular  attention  to  two  of  these 
modes,  the  sausage  mode  (m  =  0,  where  m  is  the  azimuthal  mode  number;  and 
hose  mode  (m  =  1).  These  modes  usually  have  larger  growth  rates  than  modes 
with  higher  m  numbers,  and  tend  to  cause  gross  disruption  of  the  beam.  The 
models  we  have  developed  can,  however,  be  used  to  look  at  modes  with  any 
value  of  m,  and  we  give  some  preliminary  results  for  m  =  2  and  m  =  3. 


I.  MODELS 

Our  approach  is  to  use  the  cold-fluid  equations  to  describe  the 
beam.  This  prevents  us  from  looking  at  thermal  effects  and  we  recognize 
that  present  methods  of  intense  beam  generation  tend  to  produce  beams  with 
a  significant  thermal  spread.  We  note  however  that  treatments  of  the  re¬ 
sistive  hose  instability  to  date  1-4  have  emphasized  the  effect  of  a  spread 
in  betatron  frequency  (which  is  present  in  our  cold-fluid  model)  and  that 
the  transverse  thermal  velocity  does  not  appear  explicitly  in  the  disper¬ 
sion  relations.  Thus,  the  importana  of  a  thermal  spread  as  distinct  from 
a  spread  in  betatron  frequency  is  so  far  unclear.  (We  are  presently  devel¬ 
oping  a  model  which  will  be  able  to  address  this  issue.) 

Apart  from  neglecting  thermal  effects,  our  present  model  makes  no 
additional  approximations.  It  therefore  allows  us  to  look  at  resistive 
beam  modes  in  a  consistent  manner,  without  requiring  ad  hoc  assumptions  to 
take  account  of  radial  inhomogeneity.  For  the  most  part,  we  have  used  the 
code  GRADR5;  modified  to  include  a  scalar  background  conductivity,  to  solve 
the  complete  fluid  equations.  GRADR  computes  the  linear  eigenvalues  and 
eigenfunctions  of  any  self-consistent  equilibrium  of  a  cold  single-specie 
beam.  In  addition,  we  have  made  use  of  analytical  dispersion  relations  ob¬ 
tained  from  the  fluid  equations.  These  dispersion  relations  apply  to  beams 
with  square  density  profiles  and  are  obtained  by  making  the  following 
assumptions,  commonly  employed  in  treating  resistive  instabilities. 


(v/Y)1/2  «  1  , 


(la) 


u  <<  4  iro 


(lb) 


3 


wa  «  c 


ka  «  1 


where  v  Is  the  Budker  parameter,  y  is  the  beam's  relativistic  factor,  k  is 
the  z  wavenumber  and  a  is  the  beam  radius.  The  analytical  dispersion  re¬ 
lations  have  the  following  form,  taken  from  Weinberg3, 


Jm(qa/n)\  +  r2  =  (qa)  -  am(qR  )Jm(qa)  (2) 


m  (qa/n) 


mi  2  2  1/2 

where  J  ,  H‘  '  are  Bessel  functions  of  order  m,  q  *  (4n iaaia  /c  )  with 
m  m 

Im(q)  >  0,  am(qR)  =  Hm- 1 ^CIR ^ /Jm_ i (C1R )  for  ,n  >  °»  where  R  is  the  Plasma 
channel  radius,  aQ  =  Hp^(qR)/J ^(qR)  and 


U-2-m)  (x+Z- ml 


*2(x) 


4m 

(X-m)  (2+m-x ■ 


i 

for  a  single  beam  specie.  In  these  expressions,  A  =  n/ug  =  (u-kv)/ug 
where  ug  is  the  beam  betatron  frequency  and  v  »  c  is  the  beam  z  veloc¬ 
ity.  Two  forms  of  Eq.  (2)  were  employed,  one  for  a  single  beam  specie 
using  Eqs.  (3),  and  one  for  two  counter-rotating  species  using  = 

\  Cn  (X )  +  n  (-A )],  \  U2(x)  +  C2(-A)]  in  place  of  n  and  £2. 

We  see  that  the  dispersion  relation  [Eq.  (2)]  depends  only  on  the  dimen¬ 
sionless  quantities  qa  and  A.  In  what  follows,  we  give  our  solutions  to 
Eq.  (2)  and  the  GRADR  dispersion  relation  in  terms  of  u/ug,  kc/og  and 
specify  the  dimensionless  parameter  4noa2ug/c2.  This  is  equivalent 
to  obtaining  the  solutions  in  terms  of  qa  and  A. 

The  modified  version  of  GRADR  was  checked  against  the  single 
specie  dispersion  relation  [Eqs.  (2)  and  (3)]  in  the  regime  where  condi¬ 
tions  (1)  are  satisfied.  The  frequencies  obtained  agreed  to  within  1%. 
GRADR  is,  of  course,  not  restricted  to  this  regime. 


II. 


RESULTS  FOR  HOSE  (m=l)  INSTABILITY 


A.  Step- function  current  profile:  3  models. 

We  compare  the  dispersion  curves  obtained  from  the  cold  fluid 
single-  and  two-species  models  and  the  rigid  beam  model  (see,  for  example, 
Ref.  4)  in  Fig.  1,  for  a  regime  in  which  (la)  -  (Id)  are  well  satisfied. 
(The  rigid  beam  dispersion  relation  is4 


1 


jJ(qa) 


Hp)(qa) 

Jj(qa) 


H^1)(qR/a) 
“  J0  (q'R/'a) 


All  the  models  agree  at  low  frequency  ( j  qa  J  -*•  0),  where  the  dispersion  re¬ 
lation  becomes 


k2v2  .  -  \  q2a2(*n  R/a  +i)u2 


(4) 


in  each  case.  As  one  would  expect,  the  graph  of  y  (the  imaginary  part  of 
w)  for  the  single-specie  model  is  not  syrmietrical  about  k  =  0.  The  disper¬ 
sion  relation  for  this  model  has  a  singularity  at  k  =  k*,  w  =  0,  and  this 
is  an  accumulation  point  for  an  infinity  of  modes  with  k  >  k*  (only  the 
hose  mode  exists  for  0  <  k  <  k*).  It  is  not  clear  what  happens  to  the  hose 
mode  at  k  =  k*.  However,  the  curve  shown  in  Fig.  1(b)  appears  to  have  the 
correct  continuation  of  the  mode  for  k  >  k*. 

The  single-specie  model  agrees  well  with  the  rigid  beam  model  as 
to  the  maximum  growth  rate.  Both  disagree  with  the  two-species  model  by 
about  a  factor  of  2  (cf.  Table  I).  This  Is  in  contrast  to  recent  work1 
which  predicts  that  beam  rotation  has  little  effect  on  the  hose  mode. 
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Rigid  Beam 

One  Specie 
Two  Species 


Fig.  1(a) 


Pig.  1(b) 

Figure  1.  Comparison  between  the  rigid  beam  (solid  lines),  one  specie 
(dotted  lines),  and  two  species  (dashed  lines)  models  of  the 
hose  instability  for  parameters  satisfying  (la)  -  (Id). 

The  real  and  imaginary  parts  of  the  frequency,  wr  and 
Y,  are  plotted  versus  the  z  wavenumber,  k.  There  is  an 
accumulation  point  at  k*  in  the  one  specie  model.  We 
have  taken  R/a=l  [cf.  Eq.  (2)]. 


The  difference  between  the  single-  and  two-species  result  can,  perhaps,  be 
explained  by  noting  that  the  hose  mode  is  closely  associated  with  the 
w-kv-u)g  =  0  (or  u-kv+ug  =  0)  resonance  of  the  beam.  For  the  two- 
species  case,  half  of  the  particles  are  rotating  with  frequency  -u>e  (or 
+wg)  and  respond  much  less  to  the  perturbing  fields  than  the  particles 
rotating  in  the  opposite  direction. 


B.  Effect  of  nonuniformity  in  current  profile. 

It  has  been  shown1’2  that  the  spread  in  betatron  frequency  due  to 
nonuniformity  in  the  beam  current  profile  can  significantly  affect  the 
growth  rate  and  asymptotic  behavior  of  the  hose  instability.  In  particu¬ 
lar,  while  the  rigid  beam  model  predicts  an  absolute  instability  in  the 
beam  frame,  models  with  a  spread  in  betatron  frequency1’2  predict  that  the 
instability  will  asymptotically  convect  backwards  along  the  beam.  Our 
cold-fluid  model  can,  of  course,  contain  a  spread  in  betatron  frequency,  so 
that  this  effect  should  be  present.  For  nonuniform  current  profiles,  how¬ 
ever,  we  can  only  obtain  u  as  a  function  of  k  numerically,  so  that  a  saddle 
point  analysis  of  asymptotic  behavior  would  also  have  to  be  performed  nu¬ 
merically.  We  have  decided  to  await  comparison  of  the  spectra  of  the  pre¬ 
sent  model  with  those  of  a  more  complete  model  before  attempting  to  do 
this.  For  now,  we  shall  concentrate  on  the  maximum  growth  rate,  ymax, 
obtained  from  the  linear  dispersion  relation  and  the  corresponding  group 
velocity  Vg(ymax)  with  which  the  point  of  peak  amplitude  propagates 
along. 


TABLE  1. 


Frequency  (tor,  y)  of  the  most  unstable  wavenumber  kmax,  for 

three  analytic  models. 


Rigid  Beam 

Si ngl e-Species 

Two-Species 

1.0 

1. 1 

.67 

(0.25,  0.13) 

(0.20,  0.13) 

(0.11,  0.07) 

TABLE  II. 

Effect  of  current  profile  on  frequency,  phase  velocity,  v$> ,  and  group 
velocity,  vg  of  the  most  unstable  wavenumber.  We  have  chosen 
Voa%g/c2  =  16.6,  where  o'  is  the  mean  conductivity  in  the  plasma 
channel  and  Uo  is  the  r.m.s.  betatron  frequency;  a,  and  w^a 
are  the  same  Tor  each  profile.  [c5g  =  6. 5xl0^sec"^  for  a  10  kA 
beam  with  a  =  0.5  cm  and  y  =  100.] 


Profi 1 e 

kmaxc/^e 

13 

3 

v./c 

<f> 

vg/c 

Square 

1.25 

(0.44,  0.19) 

0.35 

0.69 

Parabol ic 

1.74 

(0.66,  0.26) 

0.38 

0.72 

Gaussian 

2.2 

(0.64,  0.31) 

0.30 

0.67 

Bennett 

3.4 

(1.0,  0.44) 

0.29 

0.72 
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case  with  Gaussian  beam  current  and  conductivity  profiles.  We  varied  the 
width  of  the  latter,  while  keeping  the  peak  conductivity,  o(r=0),  con¬ 
stant.  The  results  in  Table  III  show  that  the  broadening  reduces  the  maxi¬ 
mum  growth  rate.  For  the  broadest  profile  in  this  Table  the  reduction  is 
about  25%,  and  further  broadening  has  little  effect  on  the  growth  rate. 

So  far,  we  have  looked  only  at  cases  where  the  plasma  return  cur¬ 
rent  is  negligible.  This  is  probably  a  good  approximation  for  beam  cur¬ 
rents  on  the  order  of  10  kA.  For  beam  currents  on  the  order  of  100  kA, 
however,  return  current  fractions,  f  =  -Ip/Ip  (Ip  =  plasma  return 
current,  Ip  =  beam  current)  near  80%  are  likely.  In  looking  at  the  ef¬ 
fect  of  the  return  current,  we  have  assumed  for  simplicity  that  its  profile 
is  the  same  as  that  of  the  conductivity  profile.  In  Table  IV,  we  show  the 
effect  of  changing  f,  assuming  Ip  has  the  same  shape  as  Ip.  The  pres¬ 
ence  of  a  large  return  current  does  not  have  a  very  dramatic  effect  on  the 
maximum  growth  rate.  However,  the  value  of  kmax  is  rather  more  sensitive 
to  f.  For  f  =  0,  the  instability  has  shifted  to  larger  wavenumbers  (cf. 
Fig.  2)  which  do  not  satisfy  condition  (Id)  very  well,  e.g.,  kmaxr0  = 

0.36  where  r0  is  the  width  of  the  beam  current  profile  while  for  f  =  0. 8, 
we  have  kmaxr0  =  0.13.  Condition  (Id)  is  commonly  used  1»3  in  treat¬ 
ments  of  resistive  modes  to  justify  the  neglect  of  all  but  the  z  component 
of  the  perturbed  beam  current  density. 

The  presence  of  return  current  causes  the  beam  to  be  unstable  at 
k  =  0  with  a  purely  imaginary  frequency,  so  that  the  whole  beam  displaces 
sideways  and  may  be  completely  ejected  from  the  plasma  channel.  This  ef¬ 
fect  is  due  to  the  repulsion  between  the  beam  and  plasma  currents.1 

If  the  conductivity  profile  is  broadened  so  that  some  of  the  re¬ 
turn  current  flows  outside  the  beam,  the  maximum  growth  rate  is  reduced,  as 
shown  in  Table  V.  The  largest  decrease  attained  is  about  30%.  Most  of 
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TABLE  III. 


Effect  of  broadened  conductivity  profiles  on  hose  instability, 
is  the  same  for  each  case.  As  in  Table  11,  we  have  4iroa2£Tg/c2 


d(r2)/fod(r2) 


^max0/ 

2.2 

(0.64,  0.31) 

2.1 

(0.54,  0.26) 

2.2 

(0.60,  0.25) 

2.2 

(0.58,  0.23) 

TABLE  IV. 


Effect  of  return  current  fraction  f.  Return  current  has  same  Gaussian 
profile  as  beam  current,  and  4n'aa2uo/c2  =  5. 24x1 02.  [wg  = 

2.  lxl()10sec"l  for  100  kA  beam  with  a  =  0. 5  cm  and  >  =  100.  T 


*(maxc/u6 


U)/U) 

P 

Vc 

(0.15,  0.10) 

0.07 

• 

o 

•» 

• 

o 

0. 

(0.15,  0.14) 

0.2 

(0.,  0.066) 

0. 

In  Table  II,  we  show  the  effect  of  increasing  the  spread  in  beta¬ 
tron  frequency  by  changing  the  beam  current  profile  while  keeping  the  total 
current  and  average  current  density  constant. 

The  beam  current  profiles  in  each  case  are,  for  0  <  r  <  a, 


Square: 
Parabol ic : 
Gaussian: 
Bennett : 


^b  =  ^bo  > 

Jb  -  JbpU  *  r2/a2)  , 

Jb  =  Jbg  exp(-r2/r02},  r0  =  a/2  , 
Jb  =  Jbb(*  +  p2/rb2)'2»  rb  =  a/3  » 


and  Jb  =  0  for  a  >  a.  The  square  and  Bennett  profiles  represent  extremes 
in  that  the  former  has  a  single  betatron  oscillation  frequency,  while  the 
latter  has  a  uniform  spread  in  uig2(r)  from  0  to  wg2(o).  For  these 
examples  the  conductivity  is  assumed  to  be  proportional  to  the  beam  den¬ 
sity,  and  the  mean  conductivity  is  the  same  for  each  profile.  The  mode 
considered  is  the  lowest  radial  eigenmodes,  since  this  is  likely  to  be  the 
one  least  affected  by  thermal  effects.  We  see  that  the  maximum  growth  rate 
increases,  and  moves  to  shorter  wavelength,  as  the  profile  becomes  more 
peaked  on  axis.  The  group  velocity  changes  little. 


C.  Effect  of  conductivity  and  return  current  profiles. 

Ionization  produced  by  secondary  electrons  (delta  rays)  is  expec¬ 
ted  to  produce  significant  conductivity  outside  the  beam  channel.6  Furth¬ 
ermore,  saturation  of  the  density  of  ionized  electrons  within  the  beam  may 
result  in  a  conductivity  profile  flatter  than  that  of  the  beam  current.  To 
investigate  the  effect  of  broadening  the  conductivity  profile,  we  chose  a 
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this  decrease  can  be  achieved  with  a  small  amount  of  broadening  (e.g.,  10% 
of  the  return  current  outside  the  beam).  We  also  see  from  Table  V  that 
broadening  leads  to  a  sharp  decrease  in  the  growth  rate  of  the  k  =  0  mode. 
This  is  due  to  the  enhanced  production  of  eddy  currents  outside  the  beam, 
which,  by  Lenz's  Law  tend  to  preserve  the  unperturbed  magnetic  field  con¬ 
figuration. 
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III. 


RESULTS  FOR  SAUSAGE  (m  =  0)  INSTABILITY 


In  presenting  our  results  for  the  sausage  instability,  we  shall 
again  consider  only  the  lowest  radial  eigenmode.  We  have  looked  at  the 
next  higher  mode  and  found  it  to  be  roughly  50%  more  unstable  than  the 
fundamental.  However,  the  higher  mode  should  be  more  affected  by  thermal 
effects,  and  until  we  have  had  an  opportunity  to  consider  these  effects,  we 
shall  postpone  presentation  of  the  cold-fluid  results  for  this  higher 
mode. 


As  in  the  case  of  the  hose  instability,  we  shall  consider  beams 
in  the  10  kA  (f  =  0)  and  100  kA  (f  *  0.8)  regimes,  and  examine  the  effect 
of  varying  f  and  of  broadening  the  conductivity  and  return  current  pro¬ 
files.  Gaussian  profiles  are  used  throughout.  The  results  are  summarized 
in  Tables  VI,  VII  and  VIII.  From  Table  VI,  we  see  that  broadening  the  con¬ 
ductivity  profile  reduces  the  m  =  0  growth  rate  by  no  more  than  7%  (cf.  25% 
reduction  for  m  =  1).  Conversely,  the  presence  of  a  large  return  current 
fraction  has  a  smaller  destabilizing  effect  than  for  m  =  1,  as  can  be  seen 
by  comparing  Table  VII  with  Table  IV. 

For  the  case  of  a  beam  with  f  =  0. 8,  broadening  the  return  cur¬ 
rent  and  conductivity  profiles  leads  to  a  reduction  of  about  10%  in  maximum 
growth  rate,  as  can  be  seen  from  Table  VIII.  Note  that  this  is  consistent 
with  the  results  in  Tables  VI  and  VII,  namely,  (1)  the  weak  effect  of 
broadening  o  and  (2)  a  10%  decrease  in  ymax  as  f  goes  from  0.8  to  0. 

As  we  decrease  f,  kmax  increases  [roughly  as  (i-f)l/2],  just 
as  we  found  in  the  case  of  the  m  =  1  instability  (cf.  Fig.  3).  In  the 
present  Instance,  we  have  kmaxr0  =  0.28  for  f  *  0. 8,  and  kmaxr0  = 

0.59  for  f  *  0.  Again,  this  tends  to  violate  condition  (Id),  which  may  ex¬ 
plain  the  discrepancy  between  our  results  and  those  of  Lee,7  who  finds  that 
for  f  *  0,  beams  with  Gaussian  profiles  are  stable. 


TABLE  VI. 

Effect  of  broadened  conductivity  profiles  on  the  sausage  instability  for 
Anora^jg/c^  =  16.6  [cf.  Table  III]. 


Jod(r2)/J0d(r2) 

W/ub 

"  '  '  -  '  1 - 1 

< 

-©■ 

n 

Vg/c 

1.0 

4.7 

(1.4,0.33) 

0.3 

0. 82 

0.9 

5.0 

(1.6,0.31) 

0.3 

0.86 

0.8 

5.0 

(1.6.0.31) 

0.3 

0.  86 

0.6 

5.0 

(1.6,0.31) 

— 

0.3 

0.  86 

TABLE  VII. 


Effect  of  the  magnitude  of  the  return  current  fraction  f  on  the  sausage 
instability.  Ib,Ip,  and  0  all  have  the  same  profile,  and 
4noa2we/c2  =  5. 24x1 02  [cf.  Table  IV]. 


f 

kmax0/^ 

w/uJg 

mm 

0.8 

1.6 

(0.29,  0.17) 

0.18 

0.66 

0.64 

2.3 

(0.29,  0.16) 

0.13 

0.  68 

0.25 

3.4 

(0.29,  0.15) 

0.  08 

0.69 

0.0 

4.0 

(0.29,  0.15) 

0.07 

0.7 

0.8 

k=0. 

(0.,  0.03) 

0.0 

0.0 

0.64 

0. 

(0.,  6. 3xl0-3) 

0.0 

0.0 

0.0 

0. 

• 

o 

* 

» 

o 

0.0 

0.0 

TABLE  VIII. 


Effect  of  broadening  the  return  current  profile  (along  with  the 
conductivity  profile);  f  =  0. 8  and  4n'oa2w6/c2  =  5. 24xlQ2 
throughout,  [cf.  Table  V]. 


Ip(a)/Ip 

kmaxc/u)8 

yc 

Vg/C 

1.0 

1.6 

(0.29,  0.17) 

0.18 

0.66 

0.9 

3.0 

(0.34,  0.17) 

0.11 

0.65 

0.8 

3.3 

(0.34,  0.16) 

0.10 

0.7 

0.6 

3.6 

(0.31,  0.15) 

0.09 

0.7 

1.0 

k=0. 

(0.,  0.03) 

0. 

0. 

0.9 


0. 


0. 


0. 


0. 


As  in  the  case  of  the  hose  mode,  the  beam  has  an  unstable  sausage 
mode  at  k  =  0  for  f  *  0,  whose  growth  rate  Increases  with  f  (see  Table 
VII).  However,  the  growth  rate  is  significantly  smaller  than  that  of  the  k 
*  0  hose  mode  (cf.  Table  IV),  and  is  drastically  reduced  for  broadened  con¬ 
ductivity  and  return  current  profiles  (see  Table  VIII). 
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IV. 


MODES  WITH  HIGHER  m  NUMBERS 


To  date,  modes  with  m  >  1  have  received  little  attention  in  the 
context  of  self-pinched  electron  beams.  It  Is  thought  that  these  may  be 
strongly  damped  by  the  spread  in  betatron  frequency  of  a  radially  inhomo¬ 
geneous  beam.  We  have  not  looked  in  detail  at  m  >  1  modes  with  our  present 
model  since  the  higher  the  m  number,  the  more  important  thermal  effects  are 
likely  to  be.  However,  we  have  obtained  some  preliminary  results  for  the  m 
*  2,  3  modes  using  the  two-species  version  of  Eq.  (2).  The  results  are 
presented  in  Table  IX  and  Fig.  4  along  with  results  for  m  =  0,  1  from  the 
same  dispersion  relation.  We  have  used  Weinberg's  notation3  to  denote  the 
modes,  dividing  them  Into  four  groups,  Ap, n  B^,  Cfp  and  Cfo.  Thus, 
for  example,  ^mn  denotes  a  mode  with  azimuthal  number  m  and  radial  mode 
number  n  (n  radial  zeroes).  In  this  notation,  Aqi  denotes  the  fund¬ 
amental  sausage  mode  and  Cj  denotes  the  hose  mode.  For  a  given  value  of 
m,  the  A  and  B  series  (when  they  exist)  have  an  infinite  number  of  members, 
which  at  low  frequencies  come  together  at  accumulation  points  (cf.  Fig. 

4).  We  have  only  looked  at  the  first  two  members  in  each  case.  We  see 
that  all  the  modes  for  m  =  0-3  have  comparable  maximum  growth  rates.  Also, 
the  higher  frequency  modes  (B,  D)  have  vg/c  *  1  so  that  the  point  of  max¬ 
imum  growth  rate  convects  at  near  the  beam  velocity.  We  feel  that  these 
results  call  for  a  closer  look  to  be  taken  at  the  higher  m  number  modes 
using  a  more  complete  model  for  the  beam. 
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TABLE  IX. 


Comparison  of  modes  with  m  =  0  to  3  for  square-profile,  cold,  non¬ 
rotating  beam.  In  general,  there  are  four  sets  of  modes  for  each 
m,  Aqin  (n=l  '  “)»  Bmn»  On*  °m*  although  some  of  these  are 
missing  at  low  m-values.  We  have  Wa^/c^  =  16.6  throughout. 


Mode 

^max^e 

vA/c 

4> 

Yc 

A01 

2.7 

O 

OJ 

o 

r—4 

0.35 

0.85 

A02 

4.5 

(2.7,  0.27 ) 

0.61 

0.  90 

A11 

1.2 

(0.4,  0.15) 

0.33 

0.72 

a12 

2.2 

(1.4,  0.17) 

0.63 

0.90 

B11 

4.9 

(2.0,  0.10) 

0.41 

0.95 

Bl2 

6.4 

(3.5,  0.11) 

0.55 

0.95 

Cl 

0.66 

(0.17,  0.11) 

0.26 

0.61 

b21 

6.7 

(2.8,  0.10) 

0.42 

1.0 

b22 

7.5 

(3.6,  0.10) 

0.49 

1.0 

d2 

3.0 

(1.2,  0.08) 

0.39 

0.  95 

a31 

0.58 

(1.7,  0.12) 

2.9 

0.95 

a32 

3.2 

(4.3,  0.12) 

1.35 

0.95 

b31 

=  8.0 

=  (3.2,  0.10) 

0.4 

0.95 

b32 

=  8.0 

0.36 

1.0 

C3 

difficult  to 

track;  appears  to  join  onto  an  A  mode. 

b3 

5.1 

(2.3,  0.08) 

0.45 

0.  98 

V. 


SUMMARY  AND  CONCLUSIONS 


We  have  studied  the  growth  rates  and  group  velocities  of  the  m  = 

0  and  m  =  1  resistive  beam  instabilities  using  an  exact  cold-fluid  model 
for  the  beam  and  a  scalar  background  conductivity.  In  particular,  we  have 
looked  at  the  effects  of  the  beam  current  profile,  the  width  of  the  conduc¬ 
tivity  and  return  current  profiles,  and  the  magnitude  of  the  return  current 
fraction,  f.  For  the  most  part,  we  used  Gaussian  profiles  for  the  beam  and 
return  currents  and  for  the  conductivity,  and  we  considered  beams  in  the  10 
kA  (f  =  0)  and  100  kA  (f  *  0.8),  regimes  with  y  =  100.  We  find  that  broad¬ 
ening  the  conductivity  and  return  current  profiles  leads  to  some  reduction 
in  maximum  growth  rates  (<  10%  for  m  =  0,  <  30%  for  m  =  1).  The  main  ef¬ 
fect  of  reducing  f  is  to  shift  the  instabilities  to  shorter  wavelengths. 

We  note  some  discrepancies  between  our  results  and  those  of  other  workers, 
e.g.,  the  dependence  of  m  =  1  growth  rates  on  beam  rotation,  and  the  sta¬ 
bility  of  radially  nonuniform  beams  to  the  m  =  0  mode.  Finally,  we  have 
presented  some  preliminary  results  on  modes  with  m  =  2,  3  which  show  that 
for  a  square  beam  profile,  the  growth  rates  are  comparable  to  those  for  m  = 
0  and  1.  Work  is  in  progress  on  a  more  complete  model  with  which  we  hope 
to  resolve  the  discrepancies  mentioned,  and  also  obtain  more  accurate  re¬ 
sults  in  general  than  is  possible  with  the  present  model. 


25 


REFERENCES 


1.  H.  S.  Uhm  and  M.  Lampe,  Phys.  Fluids  23,  1574  (1980);  H.  S.  Uhm  and  M. 
Lampe,  NRL  Memorandum  Report  No.  4441  (1981). 

2.  E.  P.  Lee,  Phys.  Fluids  21,  1327  (1978);  E.  J.  Lauer,  R.  J.  Briggs,  T. 
J.  Fessenden,  R.  E.  Hester,  and  E.  P.  Lee,  ibid.,  p.  1344. 

3.  S.  Weinberg,  J.  Math.  Phys.  8,  614  (1967). 

4.  E.  P.  Lee,  Lawrence  Livermore  National  Laboratory,  UC ID- 1 6268  (1973). 

5.  B.  B.  Godfrey,  Mission  Research  Corporation,  AMRC-R-223  (1980). 

6.  N.  Carron,  private  communication. 

7.  E.  P.  Lee,  Lawrence  Livermore  National  Laboratory,  UCID-18940  (1981). 


26 


